Deformation of crosslinked semiflexible polymer networks 
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Networks of filamentous proteins play a crucial role in cell mechanics. These cytoskeletal networks, 
together with various crosslinking and other associated proteins largely determine the (visco) elastic 
response of cells. In this letter we study a model system of crosslinked, stiff filaments in order to 
explore the connection between the microstructure under strain and the macroscopic response of 
cytoskeletal networks. We find two distinct regimes as a function primarily of crosslink density and 
filament rigidity: one characterized by afline deformation and one by non-affine deformation. We 
characterize the crossover between these two. 

PACS numbers: 87.16.Ka, 62.20.Dc, 82.35.Pq 
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The study of biopolymer networks and gels lies at the 
heart of the understanding of the mechanical properties 
of the cytoplasm since the mechanical rigidity of the in- 
tracellular material is largely governed by the cytoskele- 
ton, a complex network of filamentous proteins, cross 
links, and other associated proteins pj. A key player 
in this cytoskeleton is F-actin, which exhibits significant 
rigidity on the cellular scale. The material properties of 
such semi-flexible polymer networks also pose complex 
and current problems in polymer physics. Many of the 
most basic questions concerning these common and im- 
portant networks, such as how they deform under stress, 
remain unanswered. The bending rigidity of such poly- 
mers introduces a new microscopic elastic parameter that 
can have consequences for the macroscopic elastic coef- 
ficients of the bulk, semi-flexible gel. This changes the 
traditional rubber elasticity model for the bulk properties 
of gels composed of cross-linked flexible polymers. 

In this letter, we examine a simple model for 
crosslinked rods that not only allows us to quantitatively 
test the relationship between the microscopic and macro- 
scopic elastic coefficients of a randomly cross-linked net- 
work, but also sheds light on the intimately related issue 
of the spatial distribution of the network strain. Among 
the most fundamental properties of polymer networks 
is the way in which they deform under stress. Since 
the classical theories of rubber elasticity [2j it has been 
suggested that this deformation is affine, i.e., that the 
strain is uniform, as for a sheared Newtonian fluid at low 
Reynolds number. While this leads to relatively good 
agreement with experiments for rubber-like networks, 
much evidence and continuing efforts concern systematic 
deviations that point to non-affine network strains at a 
microscopic scale. 

The assumption that the deformation field is affine 
down to length scales comparable to the smallest micro- 
scopic scales in the material is a great simplification that 
allows one to construct quantitative theories relating the 
macroscopic elastic constants of a gel to the microscopic 



properties of its constituent polymers. The validity of the 
affine approximation for this class of semiflexible polymer 
materials has, however, been the subject of some debate 
0- [J 0] Whether the deformation field is affine or not 
depends, of course, on length scale; clearly at the scale 
of the entire sample, all deformations are trivially affine 
when subject to simple shear. We explore whether this 
self-averaging property of the deformation field extends 
to shorter length scales for semiflexible networks. We 
show that the degree of non-affine strain is a function of 
length scale and degree of crosslinking. Specifically, we 
find that the range of non-affine strain can extend well 
beyond the mesh size, or correlation length, the typical 
separation between filaments 0. This occurs near the 
point of rigidity percolation for the network, or for highly 
flexible filaments. This results in elastic moduli governed 
primarily by bending of filaments under non-affine strain, 
consistent with 0,0. In contrast, we find that these 
networks become increasingly affine, even down to the 
smallest scales of the network, e.g. the mesh size, at 
high cross link density, high molecular weight, or for rigid 
filaments. Here, we find that the bulk elastic moduli con- 
verge to those predicted from affine theory We also 
quantify the degree of non-affine strain and show that 
this is, indeed, dependent on the length scale. 

As we focus on the zero-frequency, or static proper- 
ties of the system we may ignore the complexities of the 
network-solvent interaction. To isolate the importance of 
semi-flexibility upon network properties, we also ignore 
the complex, nonlinear response of the individual F-actin 
filaments 0, 0, ■ We instead study in detail the depen- 
dence of the bulk shear modulus and Young's modulus of 
the material upon the cross-link density of the polymer 
gel as well as the bending and extension moduli of the 
individual filaments. We model the network filaments via 
the Hamilitonian per unit length (5s) for a filament 
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The first term takes into account the extensional defor- 
mation of the filament Sl(s) as a function of arc- length, s 
with modulus The second term determines the energy 
stored in the filament due to bending: the local tangent 
of the filament makes an angle 9(s) with respect to the x- 
axis and the bending modulus of the filament is k. Note 
that both of these terms are quadratic. We do not explore 
nonlinearities, such as buckling when compressed beyond 
the Euler instability, at the scale of individual filaments. 
While such nonlinear effects are expected at increasing 
strains and finite temperature 0] , we seek to understand 
here the fundamental properties of semi-flexible networks 
and thus address only the linearized version of the prob- 
lem. We note, however, that thermal fluctuations can re- 
sult in an effective modulus fi ~ k 2 /(kTl 3 ) for a segment 
of length I. Thus, we consider fi and k to be independent 
parameters, even though they are both determined by 
single filament elastic properties and geometry at T = 
0. Clearly the full exploration of these networks at fi- 
nite temperature presents an interesting challenge; un- 
derstanding the T = mechanical properties of these 
networks, however, is the requisite first step towards this 
more ambitious program. In addition, we will point out 
specific aspects of our results that are likely to depend 
critically on our zero-temperature assumption. After all, 
it is now understood that for some mechanical proper- 
ties of networks, the zero temperature presents a singular 
limit as in the case of rigidity percolation. 

Real F-actin networks formed from an actin monomer 
solution have a complex geometry, arising from the dy- 
namic growth and branching of filaments 0, E| . For com- 
putational efficiency, we ignore such complications and 
consider static, isotropic networks of monodisperse fila- 
ments of length L. Each filament is represented by a line 
segment deposited with random position and orientation 
to a two-dimensional rectangular shear cell. Intersec- 
tions are identified as permanent, freely-rotating cross 
links, to mimic e.g. the attachments of double-headed 
myosin molecules in an ATP-deficient solution 01 ■ The 
mean distance between cross links is l c as measured along 
a filament. Deposition continues until the desired cross 
link density L/l c has been reached. 

The network is represented by the set of mobile nodes 
{xj} consisting of all cross links and midpoints between 
cross links (the latter so as to include the dominant bend- 
ing mode). The total system energy Tt({x.i}) is then 
expressed in terms of the {x^} using a discrete version 
of JIJ. Within our linearized scheme, this 7i({xi}) is a 
high-dimensional paraboloid with a unique global mini- 
mum, corresponding to the state of mechanical equilib- 
rium at T — 0. For the initially unstressed network, this 
minimum corresponds to zero deformation. Depending 
on whether we wish to measure the shear modulus G 
or the Young's modulus Y , a shear or uniaxial strain 
7 is applied across the periodic boundaries in a Lees- 
Edwards manner. This moves the global minimum to a 



new, non-trivial position, which we numerically find us- 
ing the preconditioned conjugate gradient method. The 
stored energy per unit area can then be calculated, which 
is 7 2 /2 times G or Y || within our linear approximation, 
and hence the network modulus can be extracted. This 
procedure is repeated for different network realizations 
until a reliable estimate of G or Y has been attained. 
Further simulation details can be found in (llj . 

Apart from the system size, there are three length 
scales in the problem: two geometric lengths L and 
l c , and a third material length scale deriving from the 
stretching and bending modulii, lb = \fnf\x. On dimen- 
sional grounds, the modulii can be written in the form 
G = 4/ (L, l c , lb) (with a similar expression for Y), pro- 
vided that the system size is sufficiently large that finite 
size effects can be ignored, as applies to all of the results 
presented here. 




FIG. 1: (Color online) Superimposed examples of the en- 
ergy distribution throughout networks of cross link densities 
(left) L/lc « 8.99 and (right) L/l c « 46.77 with l b /L = 0.006 
under the given shear strain. The thickness of each segment 
is proportional to the energy density per unit length, with 
a minimum so that all filaments are visible, and the calibra- 
tion bar shows what proportion of the total energy is due to 
stretching. 

Results are given in Fig. [5] for the shear modulus G 
and the Poisson ration v = 2Y/G — 1. G monotoni- 
cally increases with the cross link density L/l c and the 
ratio lb/L. Both G and Y simultaneously vanish at the 
rigidity percolation transition [Z/Zjcrit ~ 5.9, irrespec- 
tive of lb- This is consistent with the more precise value 
[Z/Zc]crit ~ 5.932 found by Latva-Kokko et al. 0] using 
the combinatorial "pebble game" method. We find that 
G and Y scale near the transition as ~ (L/l c — [Z/Z c ]crit) 
with / = 3.0 ± 0.2, placing it in a distinct universality 
class from both central-force rigi dity and bond bending 
without free rotation at cross links 13] . The variation of v 
is more subtle if at all, varying from ks 0.5 for high densi- 
ties to ~ 0.35±0.1 near the critical point (for comparison, 
standard stability considerations require — 1 < v < 1 in 
2D). Further details of the scaling behavior near the tran- 
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sition will be presented elsewhere [ll|. We remark that 
at finite temperatures G will remain non-zero above the 
conductivity percolation transition at L/l c w 5.42, but 
it is not clear if this small-G behavior is experimentally 
observable for the macromolecules under consideration 
here. 
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FIG. 2: Dimensionless shear modulus GL//1 versus the dimen- 
sionless cross link density L/l c for different £t>, demonstrating 
convergence to the affine solution at high densities or molec- 
ular weights. The error bars are no larger than the symbols. 
(Inset) The Poisson ratio v for lb/L = 0.006. 



Away from the critical point, G decreases as the bend- 
ing modulus k decreases. Indeed, in the limit n — > 0, 
or cquivalently lh — > 0, all filaments freely bend and the 
model reduces to a random network of Hookean springs, 
which is already known to have a vanishing G for all finite 
cross link densities ^ij. In the opposite limit lu — » oo, 
bending along a filament's length becomes prohibitively 
expensive and response becomes dominated by stretching 
modes. The same is true of the limit L/l c — > oo with lb/L 
fixed, since bending modes of wavelength greater than l c 
require an area of network to twist rather than a sin- 
gle filament, and are suppressed. This transition from 
low L/l c , bending-dominated networks to high L/l c , 
stretching-dominated ones can be clearly seen Fig. ^ 

Pure, affine shear, being a combination of rotation and 
extension, induces only stretching and compression of fil- 
aments. Furthermore, networks dominated by stretch- 
ing modes must be approximately affine. This follows 
from the observation that, far above the rigidity tran- 
sition, there can be no orientational disorder without 
a corresponding cost in bending energy. It is straight- 
forward to derive analytic expressions for the modulii 
under an affine strain. A rod of length L lying at an 
angle 9 to the x-axis will undergo a relative change in 
length SL/L = 7 xy sin#cos# in response to an affine 
strain field 7 xy . According to 0} with 69 = 0, the en- 
ergy cost (after uniformly averaging over all angles 9) 
is 6H = /j,Lj^ y /16. To calculate G, we need to express 
the number of rods per unit area N as a function of L 



and l c . The exact expression is easy to derive, but for 
current purposes it is sufficient to use the approximate 
relation L/l c w (a-l)/(l-2/a) with a = 2L 2 N/tt, valid 
for L/l c w 5 and greater (note that the full expression 
is monotonic). A further correction removes the dan- 
gling ends of the rods by renormalizing the rod lengths 
to L - 2l c . Then using G = 2H/^, 

so that G ~ (%/W)(i/l c as L/l c — > oo. Expression J2J is 
plotted in Fig. Inland gives a reasonable approximation to 
the data, with the agreement improving as L/l c increases. 
The same calculations can be repeated to give Yaffmc = 
3G a ffine and hence a Poisson ratio v a m no = 0.5, which 
is also plotted in the figure. Perhaps surprisingly, we 
find remarkably good agreement between the measured v 
and the affine value, even close to the rigidity transition, 
where the affine approximation fails. 

We observe that there is an apparent crossover from 
a bending-dominated, non-affine regime for cither high 
molecular weight (L) or for high density. The natural 
measure for this is the ratio of filament length to av- 
erage distance between crosslinks. This can be under- 
stood by the fact that, unlike networks of flexible poly- 
mers, where segments along a single polymer between 
crosslinks appear to behave as effectively independent 
network strands, the segments of semificxible filaments 
between crosslinks act in series. Thus, segments far from 
free ends are forced to deform nearly affinely by the 
many constraints on their neighboring segments imposed 
by crosslinks. This suggests a physical picture in which 
non-affine deformations are primarily associated with less 
constrained free ends. Let A denote the range of such 
non-affine regions near the filament ends. We expect this 
length A (l c , Z&) to be a function of the local density of fil- 
aments (measured by l c , the distance between crosslinks) 
and the material length Then, what determines the 
degree to which the network is affine or not is the rela- 
tive size of the affine to non-affine regions, or the ratio 
L/X. We find that, sufficiently far above the transition, 
G I Gaffinc collapses to a single master curve with the em- 
pirical choice A = y/l^/hi as shown in Fi g. B) (the origin of 
this length scale is discussed elsewhere This con- 

firms the existence of a fundamental length scale A for 
non-affine deformations along the filament backbone, as 
well as two physically distinct regimes: non-affine behav- 
ior for L < A and affine behavior for L > A. In the first 
of these, we note that the modulus depends only on n 
and not on fi, signalling a bend ing dominated regime, as 
predicted by Kroy and Freyjj, |l5| . Similar results have 
been found independently [lq . Note that the crossover 
can be reached at any fixed density by varying Zb alone, 
and thus represents distinct physics from the percolation 
transition, which exists only at a specific density. Hence 
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we use absolute density, rather than its value relative to 
the transition, in our scaling function. 
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FIG. 3: The master curve of G/G a ffl nc plotted against L/X 
with A = %/li/lb- The enlarged points for L/l c « 29.09 
correspond to the same parameter values as in Fig. [I] 
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FIG. 4: Plot of the affinity measure {A6 (r)} normalised to 
the magnitude of the imposed strain 7 against distance r/L, 
for different lh/L. The value of r corresponding to the mean 
distance between crosslinks l c is also indicated, as is the solid 
line -\ {A9 2 (r)) = 1, which separates affine from non-affine 
networks to with an order of magnitude. In all cases, L/l c ~ 
29.1 and the system size was W = ^-L. 

Our findings demonstrate that the degree of affinity 
of the deformation field does depend on the length scale 
on which one looks. This is shown in Fig. ^ where the 
quantity (A9 2 (r)) is plotted against r for different values 
of L/X. (A9 2 (r)) is the square deviation of the angle of 
rotation 9 between two points separated by a distance 
r, relative to the affine equivalent. This monotonically 
decreases for increasing r, suggesting that the deforma- 
tion appears more affine when viewed on larger length 
scales. Furthermore, the deviation from affinity at the 



mesh scale r — l c is small for networks with G ~ G a ffi ne , 
and large for those with G <C Gaffine, as seen by com- 
paring Figs. 01 and Thus we can reiterate the main 
findings of our work: (1) there are two qualitatively dis- 
tinct regimes, one affine and the other non-affine; (2) the 
physics of the crossover between these is distinct from 
the rigidity transition; and (3) the crossover is governed 
by a new length scale A, where, e.g., affine behavior is 
seen for filament lengths a few times this length. Since 
this length is expected to be of order the distance be- 
tween crosslinks, real networks can be in either regime 
depending on the length distribution. 
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